function u = stable_leapfrog(u,N,mu,x_max,a,num)
if num > 0
    u1 = u;
    %u = Lax_Friedrichs_onestep(u,N,mu,x_max,a);  
    u = Lax_Wendroff_onestep(u,N,mu,x_max,a);
    for i = 2:num
        u2 = u1;
        u1 = u;
        u = stable_leapfrog_onestep(u1,u2,N,mu,x_max,a);
    end
end